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Abstract We performed the numerical simulation of quantum turbulence produced 
by thermal counterflow in superfluid 4 He by using the vortex filament model. The pio- 
neering work was made by Schwarz, which has two defects. One is neglecting non-local 
terms of the Biot-Savart integral (localized induction approximation, LIA), and the 
other is the unphysical mixing procedure in order to sustain the statistically steady 
state of turbulence. We succeeded in making the statistically steady state without the 
LIA and the mixing. This state shows the characteristic relation L — 7 2 ^ns between 
the line-length-density L and the counterflow relative velocity v n s with the quantita- 
tive agreement of the coefficient 7 with some typical observations. We compare our 
numerical results to the observation of experiment by Paoletti et a/, where thermal 
couterflow was visualized by solid hydrogen particles. 
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1 Introduction 

Quantum turbulence(QT), which consists of a tangle of quantized vortices, has been 
investigated since the thermal counterflow experiments of Vinen [1-4] half a century 
ago, while the underlying physics is far from being fully understood [5,6]. The numerical 
simulations are the useful source of knowledge about QT, because the whole dynamics 
of this system is too complicated to be described analytically. One of the powerful 
schemes of the simulation is the vortex filament model based on pioneering works by 
Schwarz [7,8]. Schwarz performed the numerical simulation of counterflow turbulence 
under the periodic boundary condition by using the localized induction approximation 
(LIA) which neglects a non-local term of the Biot-Savart integral [8]. However he could 
not obtain the statistically steady state (SSS) because the vortices lie in planes normal 
to v ns = v n — v s to prevent them from reconnecting. Therefore the unphysical mixing 
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procedure, in which half of the vortices are randomly selected to be rotated by 90° 
around the axis defined by the flow velocity, was used, and only this method enabled 
him to obtain the SSS. This failure reminds us that the LIA is unsuitable through the 
absence of the interaction between vortices. Recently, thermal counterflow in superfluid 
4 He was visualized by using solid hydrogen particles, and the velocity distribution of 
particles was observed [9]. In this experiment, two types of particles appeared. Some 
particles move straight along the normalflow with the approximately same velocity 
as normalflow. Other particles move zigzag along the direction of superfluid and has 
different velocity from superfluid. In the latter case, particles seem to be trapped in 
the core of the vortices. We try to apply our numerical results for understanding these 
observations. 

Section 2 describes the equation of motion of vortices and the method of numerical 
calculation. In Sec. 3 we show the typical numerical results of vortex tangle. Section 4 
studies the validity of the LIA by comparing the LIA calculation with the full Biot- 
Savart calculation. In the section 5 we present the velocity distribution of vortices in 
the SSS of vortex tangle and compare it with that of particles in the experimental 
observation. 



2 Equations of Motion and Numerical Simulation 

At OK the velocity so of the filament at the point s(£,t) is given by 

Hi , „. f 2(l+l-) 1 / 2 \ Hi f (si - s) x dsi M , 

s =-s / xs // ln +— / ^ ^-3— -+v s , a , 1 
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where the prime denotes the derivatives with respect to the arc length £, and k is the 
quantized circulation, ao is a cutoff parameter corresponding to a vortex core radius, Z+ 
and I- are the length of the two adjacent line elements that hold the points between. 
The first term shows the localized induction field arising from a curved line element 
acting on itself. The second term represents the non-local field obtained by carrying out 
the integral of the Biot-Savart integral along the rest of the filament. The third term 
v s ,a is an applied field. The LIA used in some works (e.g. ,[7, 8, 10, 11]) means neglecting 
the second non-local term. In contrast the calculation without the LIA is called the 
full Biot-Savart calculation. At finite temperatures, the velocity s is given by 

s = s + as x (v n - so) - as x [s x (v n - so)], (2) 

where a and a are the temperature-dependent friction coefficients, v n is the nor- 
malfluid velocity, and so is calculated from Eq.(l). In the work [8] the author neglected 
the third term of Eq.(2). We performed the full Biot-Savart calculation in this work. 
The concrete reconnection procedure used in this work is the following. Every vortex 
is represented by a string of points at intervals of almost 5£. When a point on a vortex 
approaches another point on another vortex more closely than the fixed space resolu- 
tion A£, we join these two points and reconnect the vortices [12]. For the integration 
of the motion of Eq.(2) in time we used the classical 4th order Runnge-Kutta method. 
We usually start with an initial vortex configuration of six vortex rings as shown in 
Fig. 1 (upper left). 
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3 Simulation of Counterflow Turbulence 



In this section we present the numerical simulations of counterflow turbulence at the 
temperature T=1.9K, in the computing box 0.1 x 0.1 x 0.1 cm 3 , at applied conterflow 
velocity v n s =0.286, 0.381, 0.572 cm/s. A typical result is shown in Fig.l. Here the 




initial configuration of vortex loops evolves in the periodical box to a highly chaotic 
vortex tangle. The vortex line density L(t) is defined as the vortex line length per unit 
volume. In Fig. 2 we depict the time evolution of quantity L(t). It is seen that the vortex 
tangle goes to the SSS after growth period. In Fig. 2(b) we found out that line length 
density satisfied the characteristic relation L = j^v^s which have been obtained in 
previous experiments [13]. This relation is derived from the Vinen's equation [3], and 
also obtained by Schwarz using the LIA and the dynamical scaling [8] . We obtain 7 « 
139s/cm 2 which quantitatively agrees with the experimental observation 7 « 130s/cm 2 
[14]. 



4 Validity of the LIA 

In the similar works [8,10] the authors did not obtain the SSS of turbulence (the work 
[8] created the SSS only with the mixing.) The SSS was realized in the work [11] by 
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Fig. 2 Vortex line density as a function of time for different driving normalfruid v 

equal to 0.286, 0.381, 0.572 cm/s 



using the LI A, in which the authors mentioned that failures of previous works were 
due to the unsuitable reconnection procedure. We will discuss the main reason why 
the SSS was not realized in previous works. In order to consider the validity of the 
LIA, we compare two calculations. One uses the LIA [Fig. 3 left], and the other uses 
the full Biot-Savart law [Fig. 3 right]. We run both calculations at the temperature 




Fig. 3 Side view of vortex configuration by the LIA calculation (left) and by the full Biot- 
Savart law (right) at t=34.5s. The system is a (0.2cm) 3 cube. Applied normalfluid velocity 
Vns = 0.367cm/s. 



T=1.6K, in the computing cubic box 0.2 x 0.2 x 0.2 cm 3 , and applied counterflow 
velocity v n s=0.367cm/s. To explain the difference between the results of the LIA and 
full Biot-Savart law we introduce the dimensionless anisotropy parameter [8] 

I\\ = ^ z ^[l-(s'-f||) 2 ]^. (3) 

Here fn stands for a unit vector parallel to the v ns direction, and Q is the sample 
volume. An isotropic tangle yields J|| = 2/3. At the other extreme, if the tangle consists 
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entirely of curves lying in planes normal to v ns , 1^ = 1. The time evolution of L(t) and 
I\\(t) is shown in Fig 4. Figure 3 (left) and Fig. 4 (b) show that the many vortices lie 
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Fig. 4 Comparison of a vortex line density L(t) (a) and an anisotropy parameter I\\(t) (b). 



in the planes normal to v ns , the vortex tangle being degenerate. Since the dense part 
of vortices catch other vortices moving freely, vortices become to huddle in periodical 
planes as shown in Fig. 3 (left). This ill behavior comes from the mutual friction, which 
tends to expand vortices perpendicular to the v ns , so that vortices gradually lie in 
planes normal to v ns . However, a non-local term of the Biot-Savart integral could 
yield the velocity in a direction parallel to v ns even when vortices align in a plane 
perpendicular to v ns , thus destroying the ill structure. Hence, the calculation with the 
full Biot-Savart law can sustain the SSS in contrast to that with the LI A. 



5 Velocity Distribution of Vortices 

For comparison with experimental results [9], we present the velocity distributions of 
vortices in the SSS. We obtained the distributions as shown in [Fig. 5 left] by measuring 
the z-component of the velocity s at each points on the vortex filaments from the 
results of counterflow turbulence in which applied velocity v ns is directed to z-axis. The 
velocity s does not necessarily express the particle motion on the vortex filament which 
probably occur in the experiments. Therefore we brought this effect in our calculation. 
For the sake of simplicity, suppose the trapped particles do not have a significant 
influence on the local motion of the lines, and that viscous interaction of a particle 
with the normal fluid causes the particles to move along the filament at a rate given 
by the Stokes law with the force equal to the component of the viscous force along the 
line. With these assumption, the ^-component of the velocity s p including the effect of 
motion is derived from the similar fashion of the work [15] like 

Spz = S n • Z + (v n ~ S n ■ 2) COS 2 0, (4) 

where s n is the normal component of the vortex velocity to the vortex filament, is a 
polar angle of the vortex relative to v ns . The velocity distributions of both s z and s pz 
are shown in [Fig. 5 right]. The vortices in our simulations have lower vertical velocities 
than the superfluid velocity just like the experimental observation. For the quantitative 
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Fig. 5 (a) The vertical velocity v z distribution at temperature T=1.9K, applied velocity 
Vns =0.572cm/s. The superfluid velocity v s is shown by vertical line, (b) Vertical velocity of 
vortex tangle v z , which is a peak velocity of the distribution, as a function of v 3 . The dashed 
line corresponds to v z = v s . 



agreement with observation, we need to introduce other effects such as a distortion of 
vortices by particles. 

6 Conclusions 

The obtained numerical results with the full Biot-Savart law demonstrate that the 
initially smooth vortex rings develop to a vortex tangle of the statistically steady state. 
We compared the numerical results by the full Biot-Savart calculation to that by the 
LIA calculation. The LIA calculation could not sustain the statistically steady state in 
contrast to the full Biot-Savart calculation. We compared the velocity distributions of 
vortices obtained by our numerical simulations to the observations by the visualization 
experiment of couternow by using solid hydrogen particles. Our results agree with the 
experimental observation qualitatively in that most vortices have lower velocities than 
the superfluid velocity. 
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